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Abstract 

We have developed a formalism for microscopic Brueckner-type calcula- 
tions of dense nuclear matter that includes all types of baryon-baryon inter- 
actions and allows to treat any asymmetry on the fractions of the different 
species (n,p, A, S~, S'', S+, and E^). We present results for the different 
single-particle potentials focussing on situations that can be relevant in future 
microscopic studies of beta-stable neutron star matter with strangeness. We 
find that both the hyperon-nucleon and hyperon-hyperon interactions play 
a non-negligible role in determining the chemical potentials of the different 
species. 

PACS numbers: 26.60.+C, 21.65.+f, 13.75Ev, 21.30.-x 
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I. INTRODUCTION 



The properties and composition of dense matter at supranuclear densities determine 
the static and dynamical behavior of stellar matter [|l|-|^. The study of matter at extreme 
densities and temperatures has received a renewed interest due to the possibility of attaining 
such conditions in relativistic heavy-ion collisions at GSI, and in the near future at CERN 
and Brookhaven. 

It is believed that at extremely high densities, deconfinement will take place resulting 
in a transition from hadronic to quark matter. The transition point and its characteristics 
will depend crucially on the equation of state of matter in the hadronic phase. It is well 
known that the presence of strangeness, in the form of hyperons (A,S) or mesons (K^) will 
soften the equation of state and will delay the transition. Most investigations up to date are 
made in the framework of the mean field approach, either relativistic or non-relativistic 
with effective Skyrme interactions Microscopic theories, on the other hand, aim at 
obtaining the properties of hadrons in dense matter from the bare free space interaction. 
In this sense, Brueckner theory was developed long time ago and successfully allowed to 
understand the properties of (non-strange) nuclear matter starting from interactions that 
reproduce a huge amount of A^A^ scattering observables. A first attempt to incorporate 
strangeness in the form of hyperons within Brueckner theory was made in Refs. [§,0], latter 
extended to investigations of beta-stable nuclear matter |Tl| . A missing ingredient in these 
works was the hyperon-hyperon (YY) interaction and the results of single-particle potentials 
or binding energy per baryon with a finite amount of hyperons were simply orientative. 

The recent availability of a baryon-baryon potential ||12| covering the complete 
SU(3)xSU(3) sector has allowed to incorporate the YY potential in a microscopic calcula- 
tion of dense matter with non-zero hyperon fraction [|13[. The incorporation of all possible 
baryon-baryon interactions required the solution of the G-matrix equation in coupled chan- 
nels for different strangeness sectors: A^A^ {S = 0), YN {S = —1), YY {S = —2,-3 and 



—4). The work of Ref. |T3[ concentrated mainly on isospin saturated systems, i.e., sys- 
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terns with the same fraction of particles within the same isospin and strangeness multiplet: 
T = 1/2, 5 = (neutrons and protons), T = 0, 5 = -1 (A), T = 1, 5 = -1 (S^, S+) and 
T = 1/2, 5 = —2 (H^, In this way, the comphcations associated to different Fermi seas 
for each species of the same isospin-strangeness muhiplet were avoided and the G-matrix in 
each sector was independent on the third component of isospin. 

It is well known, however, that the presence of electrons makes nuclear star matter to be 
equilibrated against the weak /5-decay reactions for neutron fractions much larger (a factor 
of 10 or more) than that for protons |lT^- [T^j . Also, the increase of negatively charged leptons 
with the baryonic density will turn into a decrease when the appearance of negatively charged 
baryons becomes energetically more favorable. This is the case of the S~ hyperon, since 
neutralizing the proton charge with S~ instead of e~ will remove two energetic neutrons 
(pS^ ^ nn) instead of one (pe^ ^ n). It is clear, therefore, that a microscopic study of 
/3-stable nuclear matter with hyperons requires the treatment of highly asymmetric matter, 
both in the non-strange sector (protons vs. neutrons) and the hyperonic one (S^ vs TP and 
S+). In the present paper we extend the study of Ref. |TB| to allow for different fractions 
of each species. We will also explore the effect of the recently available YY interaction on 
the single-particle potential of the hyperons, a crucial ingredient to determine the baryonic 
density at which the different hyperons appear. Our aim is to present a thorough analysis 
of the properties of the different baryons in dense matter, taking into account their mutual 
interactions. We will explore different baryonic densities and compositions that are relevant 
in the study of neutron stars. 



II. FORMALISM 

In this section we present the formalism to obtain, in the Brueckner-Hartree-Fock ap- 
proximation, the single-particle energies of ra, p. A, S^, S~ and H° embedded in 
an infinite system composed of different concentrations of such baryons. We first construct 
effective baryon-baryon {BB) interactions (G-matrices) starting from new realistic bare BB 
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interactions, which have become recently available for different strangeness channels 
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A. Effective BB interaction 



The effective BB interaction or G-matrix is obtained from the bare BB interaction by 
solving the corresponding Bethe-Goldstone equation, which in partial wave decomposition 
and using the quantum numbers of the relative and center-of-mass motion (ROM) reads 



The starting energy uj corresponds to the sum of non-relativistic single-particle energies 
of the interacting baryons including their rest masses. Note that we use kinetic energy 
spectrum for the intermediate BB states. The variables k, k', k" and L, L', L" denote 
relative linear momenta and orbital momenta, respectively, while K is the linear center- 
of-mass momentum. The total angular momentum, spin, isospin and isospin projections 
are denoted by J, S*, T and Mp, respectively. As usually, Q^^{k' , K;T, M^) is the angle- 
average of the Pauli operator which prevents the intermediate baryons B and B to be 
scattered to states below their respective Fermi momenta kp and kp . This angle-average 
is shown in appendix A, together with the expressions that define the Pauli operator in a 
particular (T, Mp) channel in terms of the basis of physical states. Although we keep the 
index Mp in the bare potential matrix elements they do not really have a dependence on 
the third component of isospin since we consider charge symmetric and charge independent 
interactions. Therefore, the dependence of the G-matrix on the third component of isospin 
comes exclusively from the Pauli operator, since, as can be clearly seen in appendix A, it 



{{BsBi)k"KL"S"iJ)T, Mt\ G{uj) \ {BiB2)kKLS{J)TMT) 




T 



) 



X {{BB)k'KL'S\J)TMT\G{uj) \{BiB2)kKLS{J)TMT) , 



(1) 
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acquires a dependence on Mt when different concentrations of particles belonging to the 
same isomultiplet (i.e. different values for the corresponding kp^s) are considered. 

In comparison with the pure nucleonic calculation, this problem is a little bit more 
complicated because of its coupled-channel structure. Whereas for the strangeness sectors 
and —4 there is only one particle channel (NN and EE respectively) and two possible 
isospin states (T = 0, 1), in the S = —1{S = —3) sector we are dealing with the AiV(AS) 
and T,N(T,E) channels, coupled to T = 1/2 



V 



and the EiV(E5) channel in isospin T = 3/2 



G 



In the S 
T = 



T 



-2 sector we must consider the channels AA, AS, EN and EE in isospin states 



G^AA^AA G^AA^SAT G^AA^SE ^ 
G'B.N^AA GsN-^SN GsN-^HT, 

Gee-*aa Gse-^sat G'ss-»ee 



^ 

^EN-^EN ^EN-^AE ^EN-^EIi 

Gat.^en G^ae^ae G^ae^ee 

G^ES-^SAT GeS-^AE GeE-^EE 



and T = 2 



G 



EE^SE 



In addition, each box GsiBi^BsBi has a 2 x 2 matrix sub-structure to incorporate the 
couplings between (L, S) states having the same total angular momentum J. This sub- 
matrix reads 



\l = J,S = 0\G\L = J,S = 0) {L = J,S = 0\G\L = J,S = 1)^ 
^{L= J,S =1\G\L= J,S = 0) {L = J,S = 1\G\L= J,S = 1) j 

for spin singlet-spin triplet coupling {L = J, S — L — J, S — 1) and 

^ {L^ J -1,S ^1\G\L^ J -1,S ^1) {L^ J -1,S ^1\G\L^ J + 1,S ^1)^ 
^{L = J + 1,S =1\G\L = J -1,S =1) {L = J +1,S =1\G\L = J +1,S =1) ^ 

for tensor coupling {L — J — 1, S — 1 L — J + 1, S — 1). 



B. The baryon single— particle energy in Brueckner— Hartree— Fock approximation 

In the Brueckner-Hartree-Fock approximation the single-particle potential of a baryon 
Bi which is embedded in the Fermi sea of baryons B2 is given, using the partial wave 
decomposition of the G-matrix, by 
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^ J,1,S,T,Mt 



Jo 

X kKLSTMr] G^{EB,{kB,) + Es^kB,) + niB, + m^J l^i^s; kKLSTMr) , (2) 

if both types of baryons are identical, or by 

US'\kBj = 

it^^ ^ (2J + 1) I {Tb,TbM'''^M^'''^\TMt) \' r""' k'dkfik, kB,) 

J,1,S,T,Mt ° 

X {B1B2; kKLSTMrl G\EB,{kB,) + EB,{kB,) + ms, + niB,) \B^B2; kKLSTMr) , (3) 

if they are different. The labels Sbi,Sb2 (^Bi?^B2) denote the spin (isospin) of baryons 
Bi and B2, respectively, and {TbJb^^?^^ ^t^^TMt) is the Clebsch-Gordan coefficient 
couphng to total isospin T. The variable k denotes the relative momentum of the B1B2 
pair, which is constrained by 
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with C,Bi = Mb^/Mbj^- Finally, the weight function f{k, kB^), given by 



f{k, ks,) 



i tor ft S 1+^^^ , 

for \^B,kB, - {1 + U)k\ > k^F'\ (5) 

^i-"'-[^7^---;:+5"^^^^ otherw.se, 
results from the analytical angular integration, once the angular dependence of the G-matrix 
elements is eliminated. This is done by choosing appropriate angular averages for the center- 
of-mass of the B1B2 pair and for the value of fc^j which enters in the determination of the 
starting energy. See appendix B for details. 

If the baryon Bi is embedded in the Fermi seas of several baryons Bi, B2, -B3, . . ., 
including its own Fermi sea, then its single-particle potential is given by the sum of all the 
partial contributions 

UBAk)=T.uS'\k) (6) 

( B ) 

where Ub. (k) is the potential of the baryon B^ due to the Fermi sea of baryons Bj. In 
this expression k denotes the single-particle momentum of particle Bi. The non-relativistic 
single-particle energy of baryon B is then given by 

This is precisely the single-particle energy that determines the value of the starting energy 
Lj at which the G BiB2*->B3B^-^^^T^i^ in Eq. (or (^)) should be evaluated. This implies a 
self-consistent solution of Eqs. ([^), (|^) (or (|^)) and (|^. The Fermi energy of each species is 
determined by setting k to the corresponding Fermi momentum in the above expression. 

C. Energy density and binding energy per baryon 



The total non-relativistic energy density, e, and the total binding energy per baryon, 
B/A, can be evaluated from the baryon single-particle potentials in the following way 



- = - (9) 

where p is the total barionic density. The density of a given baryon species is given by 

^3 

Pb = iri = , (10) 
where xb = Pb/ P is the fraction of baryon B, which is of course constrained by 

E^B=1. (11) 

B 

III. RESULTS 

We start this section by presenting results for the single-particle potential of each baryon 
species, as a function of the baryon momentum, for several baryonic densities and various 
nucleonic and hyperonic fractions. We have focussed on results for the Nijmegen model (e) of 
the recent parametrization |T2[, since it gives, together with model (f), the best predictions 



for hypernuclear observables |T^, apart from reproducing the YN scattering scattering data 
as well as the other models. We will restrict our calculations to matter composed of neutrons, 
protons, A's and S~'s, since these last two hyperons species are the first ones to appear Jill. 



This is also confirmed on the recent study of /3-stable neutron star matter ||18[ up tp baryonic 
density 1.2 fm~'^. 

In Fig. we show our results for non-strange nuclear matter at normal density, po = 0.17 
fm~^, and three proton fractions {xp = 0.5xn, 0.25a; at and 0), where x^v is the fraction of 
non-strange baryons, which in this case is 1. We also show the hyperon single-particle 
potentials, denoted with the label "old", obtained with the Nijmegen 1989 version of the 



YN interaction [|T^. On the right panel, corresponding to symmetric nuclear matter, we 
see that neutrons and protons have the same single-particle potential, of the order of —79 
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MeV at zero momentum. Looking at the middle and left panels we see how, as the fraction 
of protons decreases, the protons gain binding while the neutrons lose attraction. This is 
a consequence of the different behavior of the A^A^ interaction in the T = and T = 1 
channels, the T = channel being substantially more attractive. The potential of the 
proton is built from more T = than T = 1 pairs and hence becomes more attractive. The 
A single-particle potential in symmetric nuclear matter turns out to be around —38 MeV at 
k = and has a smooth parabolic behavior as a function of k. This result is larger than the 
value of —30 MeV obtained when one extrapolates to large A the s-wave A single-particle 



energy of several hypernuclei pO[. It is also much larger in magnitude than the value of 
around —24 MeV pl|-p^ which is obtained using the 1989 version of the Nijmegen YN 
potential |19[ with the standard choice for the spectrum of the intermediate YN states in 



the Bethe-Goldstone equation. The value of the E single-particle potential at A; = of —20 
MeV is somewhat more attractive than that obtained with the 1989 potential of around 
— 17 MeV. The function t/s-(^) remains pretty constant in the range of momenta explored. 
Apart from the different size, the new single-particle hyperon potentials also show a totally 
different behavior with increasing asymmetry than that observed for the potentials obtained 
with the 1989 Nijmegen YN interaction. While the old A single-particle potential turns 
to be slightly more attractive with increasing neutron fraction (i.e. going from the right 
panel to the left one), the new one becomes slightly more repulsive. The changes for the S~ 
single-particle potential are more drastic. While the 1989 interaction gives a potential 
which shows a little change with increasing neutron fraction, the new potential becomes 
strongly attractive. The value at A; = for the S~ potential changes from about —20 MeV in 
symmetric nuclear matter to —37 MeV in neutron matter. This has important consequences 
in the composition of dense matter: if hyperons feel substantially more attraction, their 
appearance in dense matter will happen at lower density. We note that our results with 



the 1989 Nijmegen interaction are consistent with those shown in ||TT|, where the same 
YN interaction is used. Some differences are found in the magnitude of the single-particle 
potentials which should be adscribed to the use of a continuum spectrum prescription in the 
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case of ITTI . 



Having established how the nucleons affect the single-particle potential of hyperons it 
is necessary to investigate the influence of a finite fraction of hyperons on the hyperons 
themselves and on the nucleons. This is visualized in Figs. ^ and |^ that show the single- 
particle potentials of the different baryons as functions of the momentum. Figure |^ shows 
results at p = 0.3 fm~'^ and a hyperon fraction xy = 0.1, which is assumed to come from only 
S~ (top panels) or split into S~ and A hyperons in a proportion 2 : 1, hence Xs- = 2xy/3 
and xa = Xy/3 (bottom panels). The panels on the right correspond to symmetric proton- 
neutron composition {xp = Xn = 0.5xn, where xn = 0.9) and the ones on the left correspond 
to a higher proportion of neutrons (xp = 0.25xAr, x„ = 0.75xAr). Starting at the upper-right 
panel we observe that the presence of S~ hyperons already breaks the symmetry between 
the proton and neutron single-particle potentials in a symmetric nucleonic composition, 
the neutrons feeling around —10 MeV more attraction. This is due to a different behavior 
between the Il~n interaction which only happens via the attractive T = 3/2 channel and 
the S~p interaction that also receives contributions from the very repulsive T = 1/2 T,N 
component. In fact, the difference bewteen the neutron and proton potentials is not as 
pronounced as we move to the lower panel on the right, where some hyperons are 
replaced by A hyperons which act identically over protons and neutrons. In the upper left 
panel, where we have increased the neutron fraction in the non-strange sector, we observe 
the typical pattern for the nucleon single-particle potentials commented in Fig. ||: the 
particle with the smallest fraction (i.e. the proton) shows more binding. However, this 
behavior is partially compensated by the presence of a sea of S~ which provides attraction 
(repulsion) to the neutron (proton) single-particle potential. We also observe that the 
feels more attraction, as a consequence of having replaced some repulsive pairs 
by attractive S^n ones. The A loses binding because the Fermi sea of neutrons is larger 
and their contribution to the A single-particle energy explores higher relative momentum 
components of the effective An interaction, which are less attractive than the small relative 
momentum ones. Finally, since the Fermi sea of hyperons is small, the differences observed 
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on the potentials by going from the top panels to the corresponding lower ones (which 
amounts to replacing hyperons by A ones) are also small. 

Similar effects are found in the results reported in Fig. ^, which have been obtained 
for a baryonic density p = 0.6 fm~^, where it is expected that nuclear matter in (3 equilib- 
rium already contains hyperons |TT|. The single-particle potential of the A hyperon is less 



attractive than that for p = 0.3 fm^^ while that of the S~ is very similar. It just gains 
somewhat more attraction when the number of neutrons increase relative to that of protons 
in going from the right panels to the left ones. As for the nucleon single-particle potentials 
we observe, also on the left panels, that the attractive S^n interaction is enhanced at these 
high densities and makes the neutron spectrum more attractive than the proton one, even 
in the asymmetric situation when one would expect the protons to be more bound. 

To assess the influence of the YY interaction we represent the separate contributions 
building the A single-particle potential in Fig. ^ and those for the S~ one in Fig. |^, for 
a baryonic density of 0.6 fm~^. The hyperon fraction of xy = 0.1 is split into fractions 
Xg- = 2xy/3 and xa = xy/3 for S~ and A hyperons, respectively. The results on the 
right hand side of Figs. ^ and |^ correspond to the symmetric nuclear case and those on 
the left to a neutron fraction three times larger than that of protons. We see that the 
contribution to the A potential from the A hyperons, represented by the dash-dotted line, is 



attractive and almost negligible, due to a weak attractive AA effective interaction and 
to the small amount of A particles present. On the contrary, the contribution from the S~ 
hyperons is larger, of the order of —10 MeV in nuclear-symmetric matter and slightly less in 
nuclear-asymmetric one, which is comparable in size with the contribution from protons and 
neutrons. This example clearly shows the important role of the YY interaction in modifying 
the properties of the A hyperon. The A acquires more attraction and its appearance in dense 
matter becomes more favorable with respect to the situation in which the YY interaction was 
neglected. The fact that the neutron (thin solid line) and proton (dotted line) contributions 
to the A single-particle potential are not the same in nuclear-symmetric matter is due to 
the E~ hyperons which make the neutrons feel more attraction and, consequently, the An 
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pairs explore the effective AA^ interaction at smaller energies, where it is less attractive. 
The different contributions to the S~ potential are shown in Fig. ^. The A hyperons (dot- 
dashed line) contribute very little due to the reduced value of their Fermi momentum. The 
contribution of the S^S^ pairs (long-dashed line) is very important, of the order of —25 
MeV in symmetric nuclear matter, and becomes crucial due to the fact that the neutron 
(thin solid line) and the proton (dotted line) contributions, which amount each one to about 
50 MeV in magnitude, almost cancel each other. In the left panel, the replacement of protons 
by neutrons, lowers the S~ single-particle potential considerably, by about 25 MeV. Again, 
neglecting the YY interactions here would have made the S~ potential about 20 — 25 MeV 
less attractive, retarding its appearance in dense matter. 

The analysis of the structure of /3-stable matter requires the knowledge of the chemical 
potential (/i^) of each baryon, defined at zero temperature as the single-particle energy of the 
Fermi momentum [Eq. (^]. In Fig. ^we show the chemical potentials as functions of density 
for different nucleon asymmetries and hyperon fractions. Note that the curves are measured 
with respect to the nucleon mass and contain, in addition to the non-relativistic Fermi 
energy, the baryon mass of each species. The top panels show the results for asymmetric 
nuclear matter (x„ = 3xp = 0.75xj\f) whereas the bottom panels stand for the symmetric 
case. On the left panels we show results for purely nucleonic matter {xy = 0), on the central 
panels we have x^- = xy = 0.1, while on the right panels xy is distributed into x\ = xy/3 
and Xs- = 2xy/3. The behavior of the chemical potentials when increasing the nucleonic 
asymmetry as well as the hyperonic fraction follows closely the trends observed in Figs. 0, |^ 
and ^ for the single-particle potential at densities p = 0.17, 0.3 and 0.6 fm~^, respectively. 
We just have to consider here that the curves in Fig. ^ also contain the kinetic energy of the 
corresponding Fermi momentum. It is interesting to comment on the high density behavior 
of the chemical potentials, since this will determine the feasibility of having hyperons in beta- 
stable neutron star matter. In symmetric nuclear matter, both the A and the S~ chemical 
potentials show, from a certain density on, an increase with increasing density which is very 
mild as compared to that assumed by phenomenological YN interactions [|^. When the 
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number of neutrons over that of protons is increased (top panels), the A chemical potential 
barely changes because of the similarity between the An and Ap interaction. However, the 
S~ hyperon acquires more binding due to the dominant H^n attractive pairs over the S~p 
repulsive ones. This will favor the appearance of S~ in dense neutron star matter, through 
the nn pT,^ conversion, when the equilibrium between chemical chemical potentials is 
achieved at both sides. Once a Fermi sea of hyperons starts to build up, however, the 
neutrons become more attractive moderating, in turn, the appearance of S~ hyperons. As 
we see, the composition of dense neutron star matter in equilibrium will result from a delicate 
interplay between the mutual influence among the different species. In fact, one needs to 
find, at each baryonic density, the particle fractions which balance the chemical potentials 
in the weak and strong reactions that transform the species among themselves. This study, 
which is beyond the scope of the present work, will be presented in a separate publication 

One of the novelties of this work is that we allow for different concentrations for the 
baryon species. Therefore, we can explicitly treat the dependence of the G-matrix on the 
third component of isospin which comes from the Pauli operator of species B,B that may 
have, even when belonging to the same isospin-strangeness multiplet, different Fermi mo- 
menta. See appendix A for more details. 

In Fig. 1^ we report the diagonal SA^ —>■ UN G- matrix elements in the ^5*0 channel, as a 
functiom of relative momentum for a density p = 0.6 fm~'^, taking xa = and x^- = Xy = 
0.1. The top panels correspond to the isospin T = 1/2 channel and the lower ones to the 
T = 3/2 one. The panels on the right are for symmetric nuclear matter, = 0.5a; AT, 

while those on the left correspond to Xp = 0.25xn- The starting energy and center of mass 
is the same for all the curves shown in the same plot, thus the dependence on Mt comes 
exclusively from the Pauli operator. Note that different pairs of particles contribute to each 
(T, Mt) combination. The case (T, Mt) = (1/2, +1/2) receives contributions from T.'^n and 
L, p pairs while J:-p and S°n contribute to {T,Mt) = (1/2,-1/2). In the case of isospin 
T = 3/2 one has contributions from E~n {Mt = —3/2), T,^p {Mt = —1/2), E+n, 
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(Mt = +1/2) and S+p (My = +3/2). We observe that the curve corresponding to the third 
component less affected by Paui blocking is always more attractive as the phase-space 
for intermediate states, which induce attractive corrections to the potential matrix elements, 
is larger. This is clearly seen in the top panel on the right, since the dotted line contains a 
channel with the S~ hyperon. When the nucleonic asymmetry is increased by going to the 
panel on the left, the effects of Pauli blocking on the neutrons are more important than those 
on the S~ hyperons. This is the reason for the full curve to appear above the dotted one, 
since the (T, My) = (1/2, +1/2) case receives contributions from TPp and S^n pairs in a 
proportion 1:2 and it contains relatively more neutrons than the case (T, Mt) = (1/2, —1/2) 
with and pairs in a proportion 2:1. In the case of T = 3/2 we observe that the 
asymmetry on the S multiplet barely induces any dependence on in the G-matrix, as 
can be seen from the bottom panel on right. However, one can observe differences when 
going to asymmetric nuclear matter on the left panel since the Pauli blocking on S~n pairs 
{Mt = —3/2) is enhanced over that on S+p pairs [Mt = +3/2). As we can see, in all cases 
considered here the dependence of the G-matrix on the third component of the isospin is 
very weak and can almost be neglected. We have also encountered this weak dependence in 
the other BB G-matrices. Therefore, a presumably good strategy and less time consuming 
would be to obtain the G-matrices in isospin saturated systems and, afterwards, calculate 
the single-particle potentials by folding the "approximate" effective interactions with the 
different baryon Fermi seas. 

We finish this section by reporting in Fig. ^ the binding energy per baryon as a function of 
density. The right and left figures describe symmetric and asymmetric {xn = Sxp = O.TSxtv) 
nuclear matter, respectively. In the top panels, we show the binding energy with x^- = for 
several values of x\ while in the bottom panels we consider x\ = and vary the concentration 
of E~ hyperons. The binding energy per baryon, calculated according to Eqs. (|),(0), is 
the result of a balance between the average kinetic energy of each baryon Fermi sea and 
the contribution from the mutual interactions, given by the average of the single-particle 
potential of each species. In order to identify the effects of the YY interaction on the binding 
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energy per baryon we have also included a curve corresponding to a calculation with a 10% 
of hyperons (either A's or S^'s) where the YY interaction is turned off (dash-dotted line). 
In both cases, turning the yF-interaction on results in a gain of binding energy which is 
larger in the case of S^. The binding energy per baryon shows a saturation density, i.e. a 
density for which the thermodynamic pressure is zero, which is too high when we consider the 
composition with only nucleons. The location of this saturation density is little affected when 
the percentage of hyperons is increased. When a small amount of nucleons is substituted 
by hyperons there is automatically a decrease of the kinetic energy contribution because the 
hyperons can be accomodated in lower momentum states and in addition have a larger bare 
mass. The analysis of the influence of the effective interaction on the binding energy must 
be made separately for A's and S~'s. Although the effective AA^ and AA interactions are 
clearly less attractive than the A^A^ one, the reduction of kinetic energy is clearly enough 
to compensate for the loss of binding energy when a 10% of nucleons is substituted by A's. 
Notice, however, that we have to consider the AA interaction in order to obtain this increase 
of binding with respect to the pure nucleonic case. At xa = 30% the loss of kinetic energy is 
not enough to compensate for the loss of attraction from the effective interactions and less 
binding energy than the case with only nucleonic degrees of freedom is obtained. Looking at 
the lower panels for the hyperons we observe that the binding energy per baryon gains 
more attraction as compared to the case for A's. This is due, essentially, to the larger loss 
of kinetic energy due to the larger mass of the E~. In general, the replacement of nucleons 
by hyperons produces a gain in binding energy and a softening of the equation of state. The 
appearance of hyperons in beta-stable matter, the softening of the equation of state and its 
implications on the properties of neutron stars are deferred to a future study |18l . 



IV. CONCLUSIONS 

In this work we have developed the formalism for microscopic Brueckner-type calculations 
of dense nuclear matter with strangeness, allowing for any concentration of the different 
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baryon species. 

By relating the Pauli operator to the different pairs of physical particles that contribute 
to the particular (T, Mr) channel (see appendix A), we have been able to obtain the Mt 
dependence of the effective interaction (G-matrix) between any two species. 

We have seen, however, that the dependence of the G-matrix on the third component 
of isospin is weak enough to allow, in future studies, for a simpler strategy consisting of 
obtaining the effective interactions in isospin saturated situations {kp'^ — kp\ kp^ ^ — 
kp — kp 1 kp = kp ). The various single-particle potentials can then be obtained 
by folding the approximate effective interactions with the Fermi seas of the different species. 

We have studied the dependence of the single-particle potentials on the nucleon and 
hyperon asymmetries, focussing on situations that can be relevant in future studies of beta- 
stable neutron star matter with strangeness. This is why, apart from neutrons and protons, 
we have only considered the A and E~ hyperons, which are the first ones expected to 
appear. We have compared the symmetric nuclear matter composition (x„ — Xp — xj^) 
with the asymmetric case containing a large fraction of neutrons {xn = 'ixp = 0.75xn), for a 
small, but relevant, hyperon fraction xy — 0.1. This fraction may be fully composed by E~ 
hyperons {x^- — xy) or contain also a small proportion of A's {xs- — 2xa — 2xy/3). We 
find that the presence of hyperons, especially S~, modifies substantially the single-particle 
potentials of the nucleons. The neutrons feel an increased attraction due to the T,~n effective 
interaction that only happens through the very attractive T — 3/2 T,N channel, while the 
protons feel a repulsion as the S~p pairs also receive contributions from the very repulsive 
1/2 EN one. 

By decomposing the A and E~ single-particle potentials in the contributions from the 
various species, we have seen the relevance of considering the YY interaction. For a baryonic 
density of 0.6 fm~^, a nuclear asymmetry of x„ = Sxp — 0.75xn and a hyperon fraction of 

= 0.1 (split into Xs- = 2xy/3 and xa — xy/3), we find that the hyperonic contribution 
to the A single-particle potential at zero momentum is of the order of —10 MeV (1/3 of the 
total C/a(0)), and that for the E" is of the order of -25 MeV (1/2 of the total [/s-(0)). 
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In the absence of hyperonic Fermi seas the A and E~ chemical potentials show a mild 
increase with increasing baryonic density. The presence of a Fermi sea of E~ hyperons slows 
down this increase, especially for the E~ chemical potential and in the case of asymmetric 
nuclear matter, due to the very attractive 7 = 3/2 "EN interaction acting on pairs. 
This will make the balance between chemical potentials in the strong nn — > E~p conversion 
easier and will favor the appearance of E~ at lower densities. 

Finally, we have studied the modifications of the binding energy per baryon in symmetric 
and asymmetric nuclear matter when some nucleons are replaced either by A or E~ hyperons. 
As expected, we observe an increase in the binding energy, which increases with density, 
mainly as a result of a decrease in kinetic energy because the hyperons can be accommodated 
in lower momentum states and have a larger mass. This effect will produce a softening in 
the equation of state that will influence the behavior of dense matter and the structure of 
neutron stars. 
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APPENDIX A: PAULI OPERATOR IN THE DIFFERENT STRANGENESS 

CHANNELS 



In this appendix we show how the Pauli operator Q^^, which prevents scattering into 
occupied BB intermediate states, acquires a dependence on the third component of isospin 
due to the different Fermi momenta of baryons B and B. The Pauli operator reads 

{1 for \aK + k\> and \3K - k\ > k^ , 
(A.l) 
otherwise 

where k and K are, respectively, the relative and total momenta of the BB pair, a = 



and (3 = . In order to solve the Bethe-Goldstone equation in partial 

mB + m-g mB + rrig 

wave representation [see Eq. (0)] we need to perform an angle average of the Pauli operator. 



which reads 



where 



, }:{cos9b + cos9~] if cos6'b + cos6'n > , 
QBBik,K) = {'' ^' ^ (A.2) 

if cos9b + cos9^ < 



1 if \aK -k\> kf^ , 

^°^^^ = ^ + (A.3) 

; otherwise 

2aKk 



and 



cos 9^ 



1 if \PK -k\> kf^ , 

P^K^ + k^-kf'" , (^-^^ 
2Wk 

Taking the following convention for the isospin states representing the particle basis 

|n) = 11/2,-1/2); b) = |l/2, +1/2) (A.5) 



|A) = |0,0) 



(A.6) 



) = |1,-1); |S°) = |1,0); 



(A.7) 
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S-=- 11/2,-1/2); =|l/2,+l/2) 



(A.8) 



it is easy to obtain the Pauli operator in the coupled-isospin basis, Q^^{k, K;T, Mt), for 
each strangeness sector. Note that in the following expressions we have only retained the 
dependence on the isospin labels. 
A.l. Strangeness 



Qnn{T = 0, Mt = 0) = - {Qpn + Qnp) 



(A.9) 



QNN{T^l,MT^-l)^Qr, 



Qnn{T = 1, Mt - 0) = - {Qpn + Qnp) 



Q^^{T^1,Mt = +1) = Q.^ 



pp 



(A.IO) 
(A.ll) 

(A. 12) 



A. 2. Strangeness —1 



Qan(T = -,Mt 



Q 



An 



(A.13) 



Qan (r=^,MT = +^) =Qap 



(A.14) 



(A.15) 



Qj:n (^T = i, Mt = = ^Qs+n + ^Q^^p 



(A.16) 



gEiv(T=|MT = -^) =Qs-n 



(A.17) 



/ 3 1\ 2 1 

QeAT f T = -, Mt = -- j = -QsOn + gQE-p 



(A.18) 
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Qetv (t = ^, = = ^Qs+n + ^Qsop (A.19) 

Qeat (t = I Mt = = Qs+p (A.20) 
A. 3. Strangeness —2 

gAA(T = 0, Mt = 0) = Qaa (A.21) 

Qsiv(r = 0, Mt = 0) = ^ (Qs-p + Qso„) (A.22) 

= 0, Mt = 0) - ^ (Qs+E- + Qeoeo + QE-E+j (A.23) 

gsiv(T = 1, Mt = -1) = Qs-n (A.24) 

gsiv(T = 1, Mt - 0) - ^ (Qs-p + Qso„) (A.25) 

Qsiv(T = 1, Mt = +1) = Qsop (A.26) 

QAE(r = 1, Mt = -1) = Qae- (A.27) 

gAE(T = 1, Mt = 0) = Qaeo (A.28) 

gAE(T = 1, Mt = +1) = Qae+ (A.29) 

Qse(T = 1, Mt = -1) = i (QsOE- + Qe-eo) (A.30) 

gsE(T = 1, Mt = 0) = ^ (Qs+E- + <5e-e+) (A.31) 

gsE(T = 1, Mt = +1) = ^ (gsOE+ + Qe+eo) (A.32) 
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gsE(T = 2, Mt = -2) = Qs-s- (A.33) 

Qss(T = 2, Mr = -1) = ^ (Qsos- + gs-so) (A.34) 

Q^^{T = 2, Mr = 0) = ^Qs+s- + ^QsoEO + ^gs-s+ (A.35) 

Qss(T = 2, Mt = +1) = ^ (gsos+ + Qs+so) (A.36) 

QsE(r = 2, Mt = +2) = gs+E+ (A.37) 
A. 4. Strangeness —3 

Qas (r = ^, Mt = = Qas- (A.38) 

Qas (t = ^, Mt = +i) = Qaso (A.39) 

gss (t = ^, Mt = = igsos- + ^Qe-so (A.40) 

Qes (t ^1,Mt^ +^) = ^Qs+s- + ^Qeoso (A.41) 

Qss (t = I Mt = = gs-s- (A.42) 

Qes (t = I Mt = -i) = ^g^os- + ^gE-so (A.43) 

gES (t = I Mr = = ^gs+s- + ^gsoso (A.44) 

gES (t = ^,Mt = +^) = gE+so (A.45) 



A. 5. Strangeness —4 
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Q^s{T = 0, Mr = 0) = ^ {Q^o^- + Qs-so) (A.46) 



Qse{T = 1, Mt = -1) = Qs-H- (A.47) 



gss(T = 1, Mt = 0) = ^ (Qsos- + <5h-so) (A.48) 



QshIT" = 1, Mt = +1) = Qhoso (A.49) 

Prom the above expressions it is easy to see that in isospin saturated matter matter (i.e. 

kp'^ = kp\ k^p ^ — k^p ^ — k^p ^ and k^f ^ = k'p' ■*) the dependence on the third component 
of isospin disappears. 
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APPENDIX B: AVERAGE OF THE CENTER OF MASS AND THE HOLE 

MOMENTA 



In this appendix we show how to compute an appropiate angular average of the center- 
of-mass momentum of the pair B1B2 and the hole momentum /e^j which enters in the 
determination of the starting energy in Eqs.(0) and (|^) . The center-of-mass momentum K 
and the relative momentum k of the pair -B1-B2 are defined in the following way: 

K = kB, + kB,, (B.l) 

— * — * 

^ = \/ = ^^^1 ~ "^^2 • (B.2) 

From the above expressions it is easy to write K and ksz in terms of the extrenal momentum 
ks^ and the relative momentum k, which is used as integration variable in Eqs. (0) and 

K = -{kB, - k) , (B.3) 
a 

kB, = -i/3kB, - k) . (B.4) 
a 

The angle average of the center-of-mass momentum is defined as 



d(cos 9) K'^{kBi,k, cos 9) 

W{kB„k) = J , (B.5) 

d{cos9) 



where K'^{kB^jk, cos 9) = + ^^ — 2^5^ cos 6^), with ^ being the angle between fc^i and 

k. The integration runs over all the angles for which Ifc^^l < k^p^\ Similarly, for the hole 
momentum we have 



d{cos 9) k'^^{k B-^, k, cos 9) 

k%,{kB,, k) = ^ J , (B.6) 

d{cos9) 



where k%^{kBi,k, cos 9) = ^(/3^A;|^ + k"^ ~ 2 jdkB^k cos 9). 

We can distinguish two cases in performing the angular integrals, < akp^^^ and 

f3kBi > ak^p^\ In the first case, we have two possibilities: < /c < ak'^p'^^ —fikB^-, for which 
all angle values are allowed, giving the result 
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(B.7) 



k^i^B^k) 



P'kl^ + k-^ 



(B.8) 



and ak^p^"* —(3kB^ < k < ak^p^^ +(3kBii which have the following upper limit in the value of 



cos^ 



2pkkB, 



cos( 



giving the result 



(B.9) 



K^{kB„k) = — 



kk + k'-^({f3kB,+kr-{ak^^^^)' 



(B.IO) 



kBo{kB^,k) = — 



P'kl^ +k'-^( {(3kB, + kf - {akf'^f 



(B.ll) 



In the second case, there is only one possibility: (3kB^—ak'^p'^^ < k < ak^p^^ +(3kBi and 
the result is the same as in the previous case for the zone akp^^^ —(3kBi < k < akp^^^ +(3kBi- 
The result for the values < k < jSkB^^—ak^p^^ is zero because is always larger than its 
Fermi sea. 

This kind of average defines an angle-independent center-of-mass momentum and a hole 
momentum (and therefore a starting energy) for each pair fc^^, so the angular integration 
in Eqs. and (|^) can be performed analytically. Nevertheless, we still require to solve 
the G-matrix equation for each pair of values fc^i and k, making the calculation much time 
consuming. In order to speed up the procedure we introduce another average, which gives 
equivalent results and saves a lot of time. For each external momentum /c^^, we will only 
need to solve the G-matrix equation for two values of the center-of-mas and hole momenta, 
which are obtained from 



/ d^kK'^{kB^,k, cost 

K^ikBj = ^ 



d^k 



(B.12) 
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j d^k k%^{kBi, k, cos( 



d-^k 



(B.13) 



by limiting the integral over the modulus of k to the two possibilities mentioned above. 
As before, we have the same cases (3kBi < akp^^^ and /^fc^i > akp^'^^ . Let's consider the 
first case. Now, when the integral over k in Eqs. ( [B.12|) and ( |B.13|) is limited to < < 
ak^p^^ —j3kB^ we have 



K'{kB,) = 



5 



(B.14) 



P^kl^ + Uak^^^^ -PkB,] 
5 



[B.15) 



whereas in the zone ak^p^^ —/SkB^ < k < ak^p^"* +(3kB^ the expressions are a little bit more 



tedious 



K\kB:, 



a 



kr'k%,+f3il + mkr' k%^ 



+ ( y + ^(4/? - 26/?^ - Q))kf^^"kl, + a'f3k^j^^W 



+ 



X 



[akf'^fpkB, + \{PkB,f - akf'\f5kB,f'^ 



(B.16) 



k'B2^kBi) — (^^^ ^^^Bi ~^ 3/3^ /c^ fc^^ 



o2,{B2)^,2 I 2,3; (52)"^; , 8/^^ ;5 

-a/5 /c), ' /c^^ + a fc^i + J^^^i 



[afcif^V/^fes, + ^{PkBj' - ak^p'^'^pkBj' 



[B.17) 



When fc^j = there exists only one zone of integration, < k < ak^p^\ and the average is 
very simple 



K^kBj = kl^ikB,) = h^p""'^" 



[B.18) 



Finally, in the second case, > ctkp^^\ there is also only one integration zone, 

PkB^—ak^p^^ < k <ak^p^^ +f3kBi, and the corresponding averages are 
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-k 



(B.19) 



(B.20) 
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FIGURES 




FIG. 1. Momentum dependence of the single-particle potentials for the different species at 
/9 = 0.17 fm~^, hyperon fraction xy = and several nucleon asymmetries. 
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FIG. 2. Momentum dependence of the single-particle potentials for the different species at 
p = 0.3 fm~^ and hyperon fraction xy = 0.1. The right panels correspond to symmetric nu- 
clear matter, Xn = Xp = O.Sxat, while the left ones are for asymmetric nuclear matter with 
Xn = Sxp = 0.75XN- In the top panels the hyperonic fraction is built exclusively from T,~ 
{xj2- = Xy) while in the bottom ones there is a fraction of A's {x\ = xy/S) and S~'s (x^- = 2xy/3) 
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FIG. 3. The same as Fig. for a baryon density p = 0.6 fm 
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FIG. 4. Separate contributions of each species to the A single-particle potential at p = 0.6 
fm~^ and hyperon fraction xy = 0.1 split into x^- = 2a;y/3 and xa = xy/3. The right panel is 
for symmetric nuclear matter (x„ = Xp = O.Sxjv) and the left one for asymmetric nuclear matter 
{xn = 3xp = 0.75xAr)- 
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FIG. 5. The same as Fig for the S single-particle potential. 
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FIG. 6. Chemical potentials of the different species as functions of total baryonic density, 
for different nucleonic asymmetries and strangeness fractions. The top panels correspond to the 
asymmetric nuclear matter case {xn = 3xp = 0.75xn), while the bottom ones correspond to 
symmetric nuclear matter (x„ = Xp = 0.5xn)- 
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FIG. 7. Diagonal 'EN G-matrix in the 5*0 partial wave as a function of the relative momentum 
at a density p = 0.6 fm~^, for the different {T,Mt) isospin channels. The right panels are for 
symmetric nuclear matter, Xn = Xp = 0.5xn, while the left ones correspond to Xn = Sxp = 0.75a; jv- 
In all cases x^- =0.1 and xa = 0. 
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FIG. 8. Binding energy per baryon as a function of the baryon density. In the top panels we set 
Xj^- = and show results for several values of x\, while the bottom panels correspond to x\ = 
and different fractions of E~'s. The panels on the right are for symmetric nuclear matter, while 
the left ones correspond to asymmetric nuclear matter (a;„ = 3xp = 0.75a; at). In the case of nuclear 
symmetric matter with 10 % of hyperons we also show a curve (dash-dotted line) where the YY 
interaction has been turned off. 
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